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Many high-temperature superconductors are highly polarizable ionic lattices where the Frohlich 
electron-phonon interaction (EPI) with longitudinal optical phonons creates an effective attraction 
of doped carriers virtually equal to their Coulomb repulsion. The general multi-polaron theory is 
given with both interactions being strong compared with the carrier kinetic energy so that the con- 
ventional BCS-Eliashberg approximation is inapplicable. The many-electron system is described by 
the polaronic t-J p Hamiltonian with reduced hopping integral, t, allowed double on-site occupancy, 
large phonon-induced antiferromagnetic exchange, J p 3> t, and a high-temperature superconducting 
state of small superlight bipolarons protected from clustering. 



PACS numbers: 71.38.-k, 72.15.Jf, 74.72.-h, 74.25.Fy 



I. INTRODUCTION 

It seems plausible that the true origin of high- 
temperature superconductivity is found in a proper com- 
bination of the finite-range Coulomb repulsion with a 
significant finite-range EPI as suggested by a growing 
number of experimental and theoretical studies In 
highly polarizable ionic lattices like cuprate superconduc- 
tors both interactions are quite strong (of the order of 1 
eV) compared with the low Fermi energy of doped carri- 
ers because of a poor screening by non- or near-adiabatic 
carriers [2|. In those conditions the BCS-Eliashberg the- 
ory Q breaks down because of the polaronic collapse of 
the electron bandwidth [H . 

The many-body theory for polarons has been devel- 
oped for extremely weak and strong EPI. In the weak- 
coupling limit this problem is reduced to the study of a 
structure factor of the uniform large polaron gas [5j . For 
strong coupling the problem is reduced to on-site [6[ or 
inter-site p], |a| small bipolarons on a lattice. A strong 
enhancement of T c was predicted in the crossover region 
from the BCS-like polaronic to BEC-like bipolaronic su- 
perconductivity due to a sharp increase of the density of 
states in a narrow polaronic band [H, which is missing 
in the so-called negative Hubbard U model. Nevertheless 
the theory of dense polaronic systems in the intermediate 
coupling regime remains highly cumbersome, in particu- 
lar, when EPI competes with strong electron correlations. 
Corresponding microscopic models with the on-site Hub- 
bard repulsion and the short-range Holstein EPI have 
been studied using powerful numerical techniques [t| [l(| ■ 

In most analytical and numerical studies mentioned 
above and many others both interactions are introduced 
as input parameters not directly related to the mate- 
rial. Quantitative calculations of the interaction matrix 
elements can be performed from pseudopotentials using 
the density functional theory (DFT) On the other 

hand, one can express the bare Coulomb repulsion and 
EPI through material parameters rather than computing 
them from first principles in many physically important 



cases [12J. In particular, for a polar coupling to longi- 
tudinal optical phonons (the Frohlich EPI), which is the 
major EPI in polar crystals, both the momentum depen- 
dence of the matrix element, M(q), and its magnitude 
are well known, |M(q)| = j(q))huio/ \/2N with a dimen- 
sionless 7(g) = y/ Aire 2 / nVLfno^q 2 , where Vt is a unit cell 
volume, N is the number of unit cells in a crystal, ujq is 
the optical phonon frequency, and K = £oo£o/( e o — e oo)- 
The high-frequency, and the static, e dielectric con- 
stants are both measurable in a parent polar insulator. 
As is well known, a two-particle bound state exists even 
in the weak-coupling regime, A < 0.5, due to a quantum 
(exchange) interaction between two large polarons form- 
ing a large bipolaron [l|(A is the familiar EPI constant 
of the BCS-Eliashberg theory). These weakly coupled 
large pairs overlap in dense systems, so that their many- 
particle ground state is a BCS-like superconductor with 
Cooper pairs (see below). 

Here the analytical multi-polaron theory is given in 
the strong-coupling regime for highly polarizable lattices 
with eo S> 1. 



II. GENERIC HAMILTONIAN AND ITS 
CANONICAL TRANSFORMATION 

The dielectric response function of strongly correlated 
electrons is apriori unknown. Hence one has to start with 
a generic Hamiltonian including unscreened Coulomb 
and Frohlich interactions operating on the same scale 
since any ad-hoc assumption on their range and relative 
magnitude might fail, 
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Here Tjj = T(m — n) is the bare hopping integral, \x 
is the chemical potential, i = m, s and j = n, s' in- 
clude both site (m, n) and spin (s, s') states, u(m, q) = 



(2N) 1 / 2 7('?) exp(iq • m), a, d q are electron and phonon 
operators, respectively, hi = c]c, is a site occupation 
operator, and H p h — X) q ^o^q^q + 1/2) is the polar 
vibration energy. 

In highly polarizable lattices with e 00 the 
familiar Lang-Firsov (LF) [l3| canonical transfor- 
mation e s is particulary instrumental with S — 
— J2q i hi [u(m, q)e? q — H.c.}. It shifts the ions to new 
equilibrium positions changing the phonon vacuum, and 
removes most of both interactions from the transformed 
Hamiltonian, H = e s He~ s , 
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where <Xy = T(m — n)XjXj is the renormalised hopping 
integral involving the multi-phonon transitions described 

with Xi = exp ^ q u(m, q)rf q — H.c. , and ft, = (J,+E p is 

the chemical potential shifted by the polaron level shift, 
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Here, the integration goes over the Brillouin zone (BZ) 
and E p — 0.647 eV in La2Cu04 0. The electron-phonon 
coupling constant is defined as A = 2E p N(0). In the 
case of 2D carriers with a constant bare density of states, 
iV(0) = ma 2 /2-kH 2 per spin, Eq.© places cuprates in the 
strong-coupling regime, A > 0.5, if the bare band mass 
m > m e (here a is the in-plane lattice constant). 



FIG. 1: (Color online) A few diagrams contributing to the 
second-order in 1/A polaron (a) and phonon (6) self-energies 
with multi-phonon vertexes. Straight and wavy lines corre- 
spond to the polaron and phonon propagators, respectively. 



ter, A(w„,k), 



III. WEAK-COUPLING REGIME 

For comparison, let us first consider the weak-coupling 
limit, where not only A < 0.5 but also the number of 
phonons dressing the carrier is small, E p /huJo <C 1. In A(cj„,k) = fegT 
this limit one can expand Xi in Eq.([5]) in powers of 7(g) 1 
keeping just single-phonon transitions so that (in the mo- 
mentum representation) 

^M(k,q) C t +qs c k , s (d q -dL q ), (4) 

q,k,s 

where £k = E(k) — p,, E(k) = — J2 m T(m) exp(im ■ k) is 
the bare band dispersion, and M(k, q) = 7(g) [E(k + q) — 
E(\t)]/V2N is the transformed EPI matrix element, 
renormalised by the Coulomb repulsion. There are no 
other interactions left in the transformed Hamiltonian 
since the bare Coulomb repulsion is nullified by the 
Frohlich EPI. 

Applying the BCS-Eliashberg formalism Q yields the 
master equation for the superconducting order parame- 



M(k, k - k') 2 D(cu n - uv)A(ov , k') 



|A(w„,,k' 



(5) 

where D(ui n — uj U ') = -/wo/[(w n — w„») 2 + ffuj 2 ] is the 
phonon propagator and u> n = nkBT(2n + 1) are the Mat- 
subara frequencies (n = 0, ±1, ±2, ±3, ...). Depending 
on a particular shape of the band dispersion, Eq.© al- 
lows for different symmetries of the order parameter since 
EPI is not local [14|. Here we confine our analysis to a 
simple estimate of T c by assuming a k-independent gap 
function, A(td n ). Then factorizing the kernel in Eq.([5j 
on the "mass shell", E(k') — E(k) — ui n > — uj n and lin- 
earizing Eq.© with respect to the gap function one ob- 
tains the familiar estimate of the critical temperature, 
k B T c « /iu;oexp[-l/(A - ^*)], where fi* = A/(l + XL) 
is the Coulomb pseudopotential. In our case the weak- 
coupling BCS superconductivity with ksT c <C huo exists 
exclusively due to the " Tolmachev-Morel- Anderson" log- 
arithm L = \n(p/haj ) > 1, if the EPI is retarded (i.e. 

hojQ < u). 
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IV. STRONG-COUPLING REGIME 

Actually the number of virtual phonons in the polaron 
cloud is large in oxides and some other polar lattices, 
Ep/hujQ ^> 1 with the characteristic (oxygen) optical 
phonon frequency Thoq < 80 meV, so that multi-phonon 
vertexes are essential in the expansion of the hopping 
operator u%y To deal with this challenging problem let 
us single out the coherent hopping in Eq.© averaging 
&ij with respect to the phonon vacuum, and consider the 



remaining terms as perturbation, H = Hq + H t 
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Here 



(0) 



describes free phonons and polarons coherently propagat- 
ing in a narrow band with the exponentially diminished 
hopping integral, tij = T(m — n) exp[— g 2 (m — n)], 



5 2 (m) = ^^ 7 (g) 2 [l-cos(q-m)], 



and 



H, 
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(8) 



is the residual polaron-multiphonon interaction, which is 
a perturbation at large A. In the diagrammatic technique 
the corresponding vertexes have any number of phonon 
lines as shown in Figll] for the second-order in H p _ p h 
polaron self-energy (E p sa —E p /2zX 2 ) and the phonon 
self-energy (S p ^ ps — xftw /z\ 2 ) [l5T |. where z is the lat- 
tice coordination number and x is the atomic density of 
carriers. Hence the perturbation expansion in 1/A is ap- 
plied if A ^> 1/ \/2z [H, Ell ■ Importantly there is no 
structural instability in the strong coupling regime since 

is ph | < m. 

The LF transformation, Eq.® is exact for any adia- 
batic ratio hioo/T(a). However, if the perturbation ex- 
pansion in 1/A is restricted by lowest orders, then it sig- 
nificantly overestimates polaron masses in the adiabatic 
regime, huj /T(a) < 1, for the case of the short-range 
(Holstein) EPI [l| (here T(a) is the nearest-neighbor bare 
hopping integral). The polaronic band narrowing factor, 
exp(— g 2 ) becomes very small for this EPI in the strong- 
coupling regime, which would eliminate any possibility of 
high temperature superconductivity and even metallicity 
of the small Hosltein polarons. 

In our case of the long-range (Frdhlich) EPI, Quantum 
Monte-Carlo simulations [17| show that the LF trans- 
formation provides numerically accurate polaron masses 
already in the zero order of the inverse-coupling expan- 
sion both in the adiabatic regime as well as in the non- 
adiabatic one for any strength of the Frohlich EPI. More- 
over, such small polarons 17 1 and small bipolarons Q are 
perfectly mobile in the relevant range of the coupling and 
the adiabatic ratio 111. 



The perturbation H p - p h has no diagonal matrix el- 
ements with respect to phonon occupation numbers. 
Hence it can be removed from the Hamiltonian in the 
first order using a second canonical transformation H = 
e S2 He~ S2 with 
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Here E n ,E n > and \n), \n ) are the energy levels and the 
eigenstates of Hq, respectively. Taking into account that 
the polaron Fermi energy is small compared with the 
phonon energy at strong coupling and / or sufficiently low 
doping [2j , one can neglect its contribution to E n i — E n w 
hwa w-q and project the second-order in 1/A Hamilto- 
nian T-L onto the phonon vacuum |0) with the following 
result 

U = -^2(tij5 aa > + p,8 i:j )c\cj- 

i,j 

^ t Knn c ms c nsC ra / s ;C n ' s ', (10) 
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ym n vti ,rp. 



ijj-i'j' i die x 



(0\[XJ(t)Xj{t) - e-^-^xlXylO), (11) 

and xj(t) is the Heisenberg multi-phonon operator ob- 
tained by replacing d q in Xj with d q exp(iuiot). Calcu- 
lating the integral, Eq. fTTj) with S — > +0 yields 
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where /(mn,m'n') = (1/2N) ]T q 7(g) 2 [cos(q-(m-n'))+ 
cos(q • (n — m')) — cos(q • (m — m')) — cos(q • (n — n'))] . 



V. POLARONIC t J p HAMILTONIAN 

All matrix elements, Eq. p2l . of the polaron-polaron 
interaction are small compared with the polaron ki- 
netic energy except the exchange interaction, J p (m — 
n) = V™ such that /(mn,m'n') = 2 3 2 (m - n). Us- 
ing Y^=iV k lk\k = -C - Info) + Ei*{y) with C « 
0.577 and Ei*(y) sa e v /y (for large y) one obtains 
a substantial J p (m) = T 2 (m) /2g 2 (m)hu>o, which is 
much larger than the polaron hopping integral, t/J p cx 
2huiog 2 e~ 92 /T(a) « 1 ia the strong coupling limit. 
Here t is the nearest-neighbor polaron hopping integrals. 
Keeping only this exchange we finally arrive with the po- 
laronic "t-J„" Hamiltonian. 



H = -}XtjjS ss > + fi5i 



hi 
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where S m = (l/2)J2s,s' c L s Tss'C m . s > is the spin 1/2 op- 
erator (t are the Pauli matrices), fi m = ^2i S hi, and 
fi = p, + X) m ^p( m ) ^ S ^ ne chemical potential further 
renormalized by H p - p h. 

There is a striking difference between this polaronic t- 
J p Hamiltonian and the familiar t-J model derived from 
the repulsive Hubbard U Hamiltonian in the limit U ^> t 
omitting the so-called three-site hoppings and EPI [l8j]. 
The latter model acts in a projected Hilbert space con- 
strained to no double occupancy. Within this standard 
t-J model the bare transfer amplitude of electrons (t) 
sets the energy scale for incoherent transport, while the 
Heisenberg interaction (J oc t 2 /U) allows for spin flips 
leading to coherent hole motion with an effective band- 
width determined by J <C t. Using the Gutzwiller-type 
approximation to remove the constraint results in an un- 
constrained t-J model also containing a band narrow- 
ing, but purely electronic rather than phononic origin 
19]. On the contrary in our polaronic t-J p Hamiltonian, 
Eq.(fT3f there is no constraint on the double on-site oc- 
cupancy since the Coulomb repulsion is negated by the 
Frohlich EPI. The polaronic hopping integral t leads to 
the coherent (bi)polaron band and the antiferromagnetic 
exchange of purely phononic origin J p bounds polarons 
into small superlight inter-site bipolarons. Last but not 
least the difference is in the " +" sign in the last term of 
Eq.(fT3"|) proportional to n m n n , which protects the ground 
superconducting state from the bipolaron clustering, in 
contrast with the "-" sign in the similar term of the stan- 
dard t-J model, where the phase separation is expected 
at sufficiently large J [20j | . 

The cancelation of the bare Coulomb repulsion by the 
Frohlich EPI is accurate up to a 1/cq correction. This 
correction produces a long-range residual repulsion of 
(bi)polarons in the transformed Hamiltonian, Eq. (TT3"l) . 
which is small as soon as eo 3> e 2 /aJ p . The residual 
repulsion results in some screening of the Coulomb in- 
teractions responsible for the doping dependence of the 
bipolaron binding energy and of the (bi)polaron mass 
In layered polar insulators the static dielectric constant 
could be anisotropic, which together with local field cor- 
rections might result in different EPI matrix elements 
for in-plane and out-of-plane polarised optical phonons, 
respectively. The difference is not considered here. 



VI. PROJECTION ONTO BIPOLARONIC 
HAMILTONIAN AND HIGH T c 

The polaronic t- J p Hamiltonian, Eq. (j!3p is analytically 
solvable in the limit of sufficiently low atomic density of 
carriers, x <C 1. Neglecting the first term in H, which 
is the polaron kinetic energy proportional to t <C J p , 
one can readily diagonalise the remaining spin-exchange 





FIG. 2: (Color online) A plane of the 3D polar lattice of an- 
ions (circles) and cations (crosses) (a) with doped carriers on 
anions bound by the polaronic exchange into four degenerate 
singlet bipolarons A,B,C and D (b). 



part of the Hamiltonian. Its ground state is an ensem- 
ble of inter-site singlet bipolarons with the binding en- 
ergy Afc = J p localised on nearest neighbor sites. Such 
small bipolarons repel each other and single polarons via 
a short-range repulsion of about J p . 

The kinetic energy operator in Eq. (|13l) connects sin- 
glet configurations in the first and higher orders with 
respect to the polaronic hopping integrals. Taking into 
account only the lowest-energy degenerate singlet config- 
urations and discarding all other configurations one can 
project the t-J p Hamiltonian onto the inter-site bipola- 
ronic Hamiltonian using the bipolaron annihilation oper- 
ators B m = 2~ 1 / 2 (c mt c m+a ! - c m ;c m+at ), where a con- 
nects nearest neighbors [7|. These operators are similar 
to the bond-order operators introduced later by Newns 
and Tsuei [2l| . which are weakly coupled in their model 
with the single-plane lattice vibrations via a nonlinear 
(two-phonon) EPI [Eq.(9) in Ref. [2lj|]. Actually it is 
well known that the nonlinear anharmonic corrections to 
EPI are small compared with the linear Frohlich inter- 
action in real three-dimensional solids, which makes the 
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two-dimensional model of Ref.[2l| of cuprate supercon- 
ductors unfeasible. 

Our strong-coupling projection is illustrated using a 
polar lattice, sketched in Fig[2K of anion-cation trian- 
gular planes (the in-plane lattice constant is a and the 
nearest-neighbor hopping distance is a/2) separated by 
the out-of-plane lattice constant c. For a zig-zag ladder- 
fragment of the lattice, FigHb, the projected bipolaronic 
Hamiltonian in the nearest-neighbor hopping approxima- 
tion is 

H b = -tJ2 B i A n + DlB n + D\A n + C\B n + Al +1 C n 

n 

+ Al +1 B n + B^_ 1 A n + C^_ 1 A n + H.c, (14) 

where A, B, C, D are annihilation operators of 
the four degenerate singlets. Fourier trans- 

formation Hi, yields four bipolaronic bands, 
E h2 (K) = -t[cos(Ka/4) ± y/l + 4sin( J ftTa/8) 4 ], 
E 3A (K) = t[cos(Ka/4) ± y/i + 4cos(ifa/8) 4 ] with the 
center-of-mass momentum hK . Expanding in powers of 
K one obtains the effective mass of these small singlets, 
to** = 10m*, where m* = 2H 2 /5t(a/2) 2 is the polaron 
mass. A similar Hamiltonian can be derived also for a 
square lattice, if next-nearest-neighbor hopping integrals 
are taking into account [22|. 

Small bipolarons are hard-core bosons with the short- 
range repulsion of the radius r = a/2 and a huge 
anisotropy of their effective mass since their inter- 
plane hopping is possible only in the second order of t 
[23| . The occurrence of superconductivity in bipolaronic 
strong coupling systems is not controlled by the pairing 
strength, but by the phase coherence among the electron 
pairs [6]. While in two dimensions Bose condensation 
does not occur in either the ideal or the interacting sys- 
tem, there is a phase transition to a superfluid state at 



in the dilute Bose gas [2J, [25| (here rib is the boson 
density per unit area). Using Eqs.([31 [7]) we obtain 
E p « 0AE c and g 2 « QA8E p /tiwo, allowing for a quan- 
titative estimate of T c (here E c — 2e 2 c/irKa 2 ). With 
typical values of a = 0.4 nm, c = 1.2nm, k = 5, the bare 
band mass to = m e , Hujq = 80 meV and the moderate 
atomic density of polarons, x = 0.1 (avoiding an over- 
lap of bipolarons) one obtains E p 0.55 eV, g 2 « 1.24, 
and T c w 205K. Importantly, the projection procedure 
of reducing Eg. (Tip) to Eq. (fT4")) is well justified since the 
ratio t/Jp m 0.1 is small and /cbT c <C J p , so that only 
the lowest singlet configurations can be included while 
discarding the others. 

In conclusion, it seems very likely that a peculiar 
cancelation of the long-range Coulomb repulsion by the 
long-range Frohlich EPI can help much in producing 
high-temperature superconductivity in doped polar in- 
sulators such as cuprates and other oxides, for instance 
BaKBiO. The polaronic t-J p Hamiltonian, Eq. (fl3|) de- 
rived here from the bare long-range Coulomb interactions 
could provide a novel avenue for analytical and computa- 
tional studies of superconductivity in complex ionic lat- 
tices since the repulsive Hubbard U model and its strong- 
coupling t-J projection do not explain high T c [26| . 



2mn J h 2 



fegTO** ln[ln(l/rif,r 2 ) 



(15) 
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